Evidence of a magnetic Bose glass in IPA-Cu(Clo.95Br .o5)3 from neutron diffraction 
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We report the single crystal study of the bulk magnetization and neutron scattering measure- 
ments on a quantum S=l/2 spin ladders system IPA-Cu(Clo.9sBro.o5)3 with quenched disorder. In 
zero field, the disordered spin liquid phase is preserved as in pure IPA-CUCI3. Due to the bond 
randomness, a different Bose glass phase was directly observed in H C <H<H' , which separates the 
spin liquid phase from the unconventional Bose-Einstein condensation phase. The observed finite 
value of boson compressibility (dM/dH) and lack of field-induced three-dimensional long range order 
are consistent with the theoretical prediction. 

PACS numbers: 75.10. Jm 05.30.Jp 71.27.+a 



A "Bose glass" (BG) is an exotic state of matter that 
emerges in systems of interacting bosons in the pres- 
ence of quenched disorder. At sufficiently low tempera- 
tures, disorder- free bosons are subject to so-called Bose- 
Einstein condensation (BEC). BEC can involve atoms 
in liquid 4 He( 1 - laser-cooled ions in magnetic traps, 2 
Cooper pairs in superconductors^ or magnons in mag- 
netic systems. 4 Due to peculiarities of Bose statistics, 
particles lose their individuality and occupy a unique 
quantum-mechanical state. The wave function of this 
condensate establishes long-range quantum phase coher- 
ence across a macroscopic sample. For repulsive bosons, 
quenched disorder disrupts the condensate and interferes 
with phase coherence. The result is a peculiar glassy 
state with only short-range phase correlations. 5,6 While 
some experimental evidence of this was found in ultracold 
atoms, 7 high-temperature superconductors,— and quan- 
tum magnets^^ none of the studies were direct. The 
key characteristic, namely the wave function of the con- 
densate disrupted by disorder on the microscopic scale, 
remained inaccessible. In this paper, we report a direct 
neutron diffraction observation of short range correla- 
tions of the BEC order parameter in a magnetic BG. 
This phase is realized in the quantum spin ladder com- 
pound IPA-Cu(Clo.95Bro.o5)3, where disorder is induced 
by random chemical substitution^ 

The disorder-free parent compound IPA-C11CI3 is a 
prototypical S = 1/2 AF spin ladder material with 
the ladders running along the a axis of the crystah 12 ' 13 
Nearest-neighbor spin interactions along the legs of each 
ladder are AF. Nearest-neighbor inter-leg correlations 
are ferromagnetic (FM). However, inter- leg coupling is 
dominated by next-nearest-neighbor interactions. These 
are formed by Cu-Cl-Cl-Cu superexchange pathways, 
and are robustly AF. As discussed in Ref. 13, magnetic 
anisotropy in this material is negligible. Zero-point quan- 
tum spin fluctuations in such Heisenberg ladder struc- 
tures destroy conventional magnetic order. The result 
is a non-magnetic "spin liquid" state. The lowest-energy 



excitations are a triplet of long-lived S — I quasiparticles 
with a minimum excitation energy A. For IPA-CUCI3, 
this energy gap is A = 1.17 meV 4^ The quasiparticles 
obey Bose statistics, and are mutually repulsive at short 
distances. Since the energy cost of creating each quasi- 
particle is at least A, the ground state is a vacuum of 
bosons. The vacuum persists in modest applied mag- 
netic fields. However, due to Zeeman effect, the gap in 
the S z = +1 magnon decreases linearly with increas- 
ing field H, and reaches zero at H c = A/ (9/13)- For 
IPA-CuCl 3 H c = 9.7 TM Once H > H c , the quasipar- 
ticle energy becomes negative, and macroscopic number 
of them are incorporated in the ground state. Since each 
carries a spin projection S z — +1, their density is equal 
to the uniform magnetization: (p) = m = {S z }. Simulta- 
neously, the emerging bosons undergo magnon BEG— ^ 
The signature of this quantum phase transition is the ap- 
pearance of spontaneous long-range staggered (AF) mag- 
netic order of spin components perpendicular to the di- 
rection of applied field. This transverse magnetization, 
written in complex form ^ = (S x ) +i{S v ), is the effective 
wave function of the Bose condensate. In IPA-CUCI3 it 
was previously directly probed by means of magnetic neu- 
tron diffraction, where the measured scattering intensity 
is proportional to the Fourier transform of the spin corre- 
lation function ; 13 ' 16 The BEC phase is characterized by a 
new set of magnetic Bragg peaks with half-integer Miller 
indexes. Their intensity is proportional to the square of 
the BEC order parameter. Note that in other experimen- 
tal realizations of BEC of magnons, such as that in thin 
films^ 7 - or 3 He 18 , the condensate wave function remains 
experimentally inaccessible. 

The best way to introduce quenched disorder in a mag- 
netic system is by chemical substitution. Most previous 
studies targeted the magnetic ions, randomly substitut- 
ing them by non-magnetic or different spin impurities. 
The problem with this approach is that it qualitatively 
alters the nature of the spin liquid state. Upon substi- 
tution, local S = 1/2 degrees of freedom are liberated in 
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FIG. 1: (color online). The excitation spectrum (false color 
plot) measured in IPA-Cu(Clo.9sBro.o5)3 at T = 1.5 K in zero 
applied field reveals sharp dispersive quasiparticles with an 
energy gap A = 1.24(1) meV. A typical energy scan (left 
panel, symbols) shows a well-defined peak with an energy 
width entirely due to experimental resolution (solid line, cal- 
culation) . 



direct proportion to the impurity concentration.— These 
free spins are the dominant contribution to bulk mag- 
netization, give rise to a divergent magnetic susceptibil- 
ity, and enable conventional long-range ordering at low 
temperatures. 20 In addition, each impurity becomes a po- 
tent scattering center. This causes a collision damping 
of magnons^i making the very quasiparticles that are 
supposed to condense poorly denned. Thus, any physics 
related to the BEC transition is wiped out. With this in 
mind, in the present study we adopted a radically differ- 
ent "softer" method. Quenched disorder was introduced 
in IPA-CuCi3 by a partial substitution of non-magnetic 
Br~ for the likewise non-magnetic Cl~. 10 This modifica- 
tion does not directly involve the spin-carrying Cu 2+ ions 
and does not add additional anisotropy term. Instead 
it affects the bond angles in the Cu-halogen-halogen-Cu 
superexchange pathways. The strength of magnetic in- 
teractions on the affected bonds is thereby modified, but 
their AF character is preserved. 

IPA-Cu(Clo.95Br .o5)3 single crystals were grown in so- 
lution as described in Ref. 11. The microscopic homoge- 
niety of Br-distribution was confirmed by single crystal 
X-ray diffraction studies, that also confirmed the Br con- 
tent to be within 10% of nominal. Bulk magnetization 
data were taken at T—hQ{) mK. The magnetic field was 
applied along the c axis of the crystal.™ Neutron experi- 
ments were carried out in the (h, k, 0) scattering plane on 
an assembly of five fully deuterated crystals with a total 
mass of 0.8 g and a mosaic spread of 0.7°. Inelastic mea- 
surements were performed on IN22 3-axis spectrometer 
at ILL, using a standard He-flow cryostat with a Be neu- 
tron filter after the sample and fixed final neutron energy 
Ef—5 meV. High- field diffraction data were collected us- 
ing a vertical field cryomagnet and a dilution refrigerator 
insert on the IN22, with a Pyrolytic Graphite (PG) filter 
after the sample and Ef=\A.l meV. The experimental 
resolution was calculated in the Popovici approximation. 



Previous bulk studies have shown that in IPA- 
Cu(Cli_ a; Br a ;)3 the spin liquid ground state remains in- 
tact up to about x — 13% Br-content. 11 In the present 
work we shall focus on the x — 0.05 material. While mag- 
netic susceptibility experiments revealed some residual 
free spins in this system, their concentration is negligibly 
smalliii From the bulk magnetization, where the para- 
magnetic contribution at T = 500 mK remains smaller 
than 0.003 per formula up to H = 8 T, we esti- 
mate that there are fewer than one free 5=1/2 spin 
for every 10 Br substitutions^ These are likely due to 
crystallographic defects: IPA-Cu(Clo. 956^.05)3 crystals 
grown from solution are systematically smaller and of in- 
ferior quality compared to those of the pure compound. 
The singlet ground state in IPA-Cu(Clo. 956^.05)3 was 
confirmed in recent mu-SR studies, where no long-range 
magnetic order was observed down to at least T = 
330 mK. 22 However, the clearest evidence of that this 
material is an immaculate spin liquid is provided by 
our recent inelastic neutron scattering experiments. The 
excitation spectrum shown in Fig. Q] was measured at 
T = 1.5 K and reveals well-defined bosonic quasiparticles 
with a spin gap A = 1.24(1) meV. Just like in the pure 
IPA-CuCl^ and PHCC^, the quasiparticle spectrum 
terminates at a critical wave vector h c ~ —0.2. Apart 
from the slightly larger gap energy, these excitations are 
almost indistinguishable from those in pure IPA-CUCI3. 
They are also at least as long-lived: knowing our energy 
resolution, we can estimate the intrinsic quasiparticle en- 
ergy width to be r < 0.03 meV. Thus, despite the struc- 
tural disorder, IPA-Cu(Clo. 956^.05)3 indeed remains a 
true quantum spin liquid. 

In external magnetic fields exceeding H c ~ 10 T, 
IPA-Cu(Clo.95Br .o5)3 becomes magnetized (Fig. [U[a), 
circles). 10 The magnetization derivative does not jump 
abruptly, as in IPA-CUCI3 under similar conditions 
(Fig. H[a), solid curves). Instead, it gradually increases 
between H c and H' ~ 11 T (Fig. [Ha), squares), and is 
roughly constant at higher fields. H c was determined by 
linear interpolations of the dM /dH curve, as shown in 
Fig. [5J H 1 was determined in power-law fits to the tem- 
perature dependence of the (0.5,1.5,0) magnetic Bragg 
peak measured upon cooling. The net magnetization and 
its derivative are to be interpreted as the density and 
compressibility of S z = +1 quasiparticles now present in 
the ground state. Our key result is that, unlike in the par- 
ent compound, in the doped material these quasiparticles 
initially fail to form a condensate, the latter only setting 
in at higher fields. In the pure system, sharp Bragg reflec- 
tions corresponding to magnetic long-range order appear 
simultaneously with bulk magnetization (Fig. [Hb) , solid 
curve) i 13 i 16 In contrast, in IPA-Cu(Clo. 956^.05)3 only 
broad peaks are observed around the propagation vec- 
tor (1/2,1/2,0). At T = 600 mK they are totally ab- 
sent below H c , appear somewhere between H c and H 
but remain barely detectable, and grow rapidly beyond 
H' (Fig. [2fb), circles). They persist to the maximum 
attainable experimental field of 13 T. Similar behavior 
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FIG. 2: (color online). Signatures of various quantum phases 
in IPA-Cu(Clo.95Bro.o5)3. The uniform magnetization (a, cir- 
cles, from Ref. 1(? ) represents the density of bosons in the 
ground state. The rate of increase of uniform magnetization 
is the boson compressibility (a, squares). The thin solid lines 
correspond to the disorder- free parent compound IPA-CUCI3. 
(b) The magnetic Bragg intensity measured at the maximum 
represents the effective BEC order parameter. The measured 
temperature dependence of H' is shown in the inset. The 
solid curve is a guide for the eye. 



is observed at other temperatures. The measured In- 
dependence of H' is plotted in the inset of Fig. 2. 

The main panel of Fig. [3] shows scans across one repre- 
sentative reflection. Here the dashed lines are the exper- 
imental resolution. The finite intrinsic peak widths un- 
ambiguously indicate that ordering is short-range. The 
corresponding correlation lengths are history-dependent. 
This is because in order to enter the BEC phase the sys- 
tem has to cross through the disordered BG state. The 
resulting pinning of magnetic domain walls or vortexes by 
the random potential makes true long-range order kinet- 
ically inaccessible. In a sample zero-field-cooled (ZFC) 
to 30 mK, at H = 13 T, the correlation lengths are ( a — 
70(6) lattice units along the ladder direction (et-axis of 
the crystal) and only ^ = 11.7(0.1) lattice units perpen- 
dicular to it (6-axis). If a sample is cooled in a 13 T field, 
the peaks are much sharper, though still broader than 
resolution, with ( a — 143(8) and ^ — 52(2). The peak in- 
tensity is correspondingly higher in the field-cooled (FC) 
sample, and shows hysteresis when the field is repeatedly 
decreased below H c , then again increased (Fig. HtK)). 
Similar hysteresis was observed in IPA-Cu(Clo. 956^.05)3 
at all temperatures between 30 and 750 mK, but never in 
pure IPA-CUCI3J 13 ' 16 History-dependent behavior and a 
difference between FC and ZFC samples are a signature 
of a magnetic glass. 
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FIG. 3: (color online). Neutron scans across the 
(—0.5,-1.5,0) point measured in zero-field-cooled (ZFC) 
and field-cooled (FC) IPA-Cu(Clo.9sBro.o5)3 samples at T — 
30 mK in H = 13 T applied field (symbols), in the BEC phase. 
The magnetic peaks are broader than experimental resolution 
(dashed lines), indicating that only short-range correlations 
of the BEC order parameter are present. The solid lines are 
Voigt function fits to the data. 



Based on bulk magnetization data, it was previously 
suggested that a BG emerges in IPA-Cu(Clo. 956^.05)3 
in magnetic fields exceeding H c ^ ~ 40 T, while for 
H c < H < H C 3 the system is in a BEC phased The 
present direct measurement of the effective BEC wave 
function correlations unambiguously show that the glass 
forms already at H c , as soon as the bosons precipitate 
in the ground state. This behavior can be qualitatively 
understood in the framework of recent theories ! 24 ' 25 The 
ground state of weakly interacting spin ladders are a lin- 
ear combination of valence-bond states composed of local 
singlet spin pairs. Due to the spatially randomized in- 
teraction strength, certain singlets may have a reduced 
gap energy. When the magnetic field is increased past 
H c , these are the first to become magnetized. Each such 
broken singlet corresponds to a S z = +1 boson localized 
by the quenched disorder. The phenomenon is similar 
to Anderson localization of fermions in a random po- 
tential. The role of the Pauli exclusion principle that 
enables the latter is played by the strong short-range re- 
pulsion between spin ladder excitations. The result is 
the magnetic BG with a finite boson density, but only 
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short-range frozen transverse spin (effective wave func- 
tion) correlations^ A lack of long-range magnetic order 
and a non-zero yet finite compressibility (dM/dH) are 
the two key feature of a BG. The latter distinguish it 
from other disordered phases of quantum spin systems: 
the incompressible Mott glasses and the random single 
phase that has divergent compressibility 2 ?". The quasi- 
1D character of IPA-Cu(Clo. 956^.05)3 undoubtedly fa- 
vors boson localization and helps stabilize the BG state. 
Indeed, only in one dimension does the BG phase appear 
for arbitrary weak disorder £ In higher dimensions the 
magnitude of disorder needs to exceed some threshold 
value to disrupt the condensate^ 

At a certain higher field, that we identify with H', 
the gaps associated with even the strongest bonds in the 
system are overcome by Zeeman energy. Beyond this 
point one expects to recover the coherent BEC phase,— 
in which the compressibility is constant and AF spin cor- 
relations are long-range. 4 Qualitatively, the field range 
H' — H corresponds to the energy difference between the 
strongest and weakest AF bonds. The same energy scale 
determines the doping-induced shift of the gap at H = 0. 
For IPA-Cu(Clo.95Br .o5)3 the latter is of the order of 0.1 
meV, and corresponds to a field range of 1 T, consistent 
with the measured H 1 — H c . 

In our experiments, beyond H' , the compressibility 
indeed levels off and AF correlations rapidly build up. 
However, the latter remain short-range. In this sys- 
tem, there is not necessarily a contradiction. We sug- 
gest that while at higher fields the BEC phase may be 
the true ground state, at low temperatures it remains 
kinematically inaccessible. As the external field is in- 
creased beyond H c , local correlated regions grow in size 
around each broken singlet that act as nucleation cen- 
ters. By the time H' is reached, the macroscopic sam- 
ple is a mosaic of uncorrelated finite size AF domains or 
a textured pattern of magnetic vortexes. At low tem- 



peratures the domain walls and/or vortexes are pinned 
by the random potential. The system is frozen in this 
short-range correlated state. In a field-cooled sample, the 
BG-BEC boundary is crossed at a higher temperature, 
where the domain walls and vortexes are more mobile. 
The result is fewer pinned defects, longer-range correla- 
tions and and sharper diffraction peaks. This interpre- 
tation allows us to reconcile the conclusions of previous 
bulk measurements 10 with the present neutron scattering 
study. For < H < H c IPA-Cu(Clo.9 S Br .o5)3 is a true 
spin liquid. For H c < H < H' the system becomes a 
magnetic BG. Long-range order in the form of BEC may 
be the ground state for H' < H < H c3 , but the system 
stays frozen in a metastable glassy state. Between H C 3 
and the saturation field H% ~ 60 T a second BG phase 
is realized, as discussed in Ref. 10. It is related to the 
saturation transition, that can also mapped on BEC. 28 

In summary, the unique ability of neutron scattering 
to probe the effective condensate wave function correla- 
tions and of magnetization measurements to probe boson 
density and compressibility in magnetized quantum spin 
liquids, allowed us to directly observe the exotic magnon 
BG phase. 
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